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ABSTRACT 

We present an analysis comparing the properties of 45 giant molecular clouds (GMCs) in M33 and 
the atomic hydrogen (H I) with which they are associated. High-resolution VLA observations are 
used to measure the properties of H I in the vicinity of GMCs and in regions where GMCs have not 
been detected. The majority of molecular clouds coincide with a local peak in the surface density of 
atomic gas, though 7% of GMCs in the sample are not associated with high-surface density atomic 
gas. The mean H I surface density in the vicinity of GMCs is 10 il/0 pc"^ and tends to increase 
with CMC mass as Shi oc Mq^'q. Thirty- nine of the 45 H I regions surrounding GMCs have linear 
velocity gradients of ^ 0.05 km pc~^. If the linear gradients previously observed in the GMCs 
result from rotation, 53% are counterrotating with respect to the local H I. And if the linear gradients 
in these local H I regions are also from rotation, 62% are counterrotating with respect to the galaxy. 
If magnetic braking reduced the angular momentum of GMCs early in their evolution, the angular 
velocity of GMCs would be roughly one order of magnitude lower than what is observed. Based on our 
observations, we consider the possibility that GMCs may not be rotating. Atomic gas not associated 
with GMCs has gradients closer to 0.03 km s~^ pc~^, suggesting that events occur during the course 
of CMC evolution that may increase the shear in the atomic gas. 
Subject headings: galaxies: individual (M33) — galaxies: ISM — ISM: clouds 
dynamics — ISM: molecules — radio lines: galaxies 
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1. INTRODUCTION 

Observations of molecular clouds at millimeter wave- 
lengths show that they often have systematic velocity 
gradients. If, as many authors have argued, these veloc- 
ity gradients are indicative of cloud rotation (e.g., Kutner 
et al. 1977; Blitz 1993; Rosolowsky et al. 2003), then an- 
gular momentum conservation should contribute to our 
understanding of the formation and evolution of molec- 
ular clouds. Significant progress was made in this direc- 
tion when Rosolowsky et al. (2003) studied the angular 
momentum properties of giant molecular clouds (GMCs) 
in M33. Their measurements, based on high-resolution 
^^CO(J = 0^-1) observations, showed that M33 molec- 
ular clouds have velocity gradient magnitudes that are 
comparable to Milky Way GMCs (~ 0.1 km s~^ pc~"^). 
They found that typical gradients of the GMCs are 5 
to 10 times smaller than would be expected from simple 
formation theories. Furthermore, just as Galactic GMCs 
are sometimes observed with gradients that are not par- 
allel to the sense of Galactic rotation (e.g.. Blitz 1993; 
Phillips 1999; Imara & Blitz 2010), Rosolowsky et al. 
(2003) found that the gradient directions of M33 GMCs 
are often not in alignment with the galaxy rotation. This 
is the so-called "angular momentum problem" of GMC 
formation. 

It is our goal in this paper to address the question 
of the origin of velocity gradients in GMCs by doing a 
detailed analysis of the gradients in the atomic gas as- 
sociated with the M33 molecular clouds catalogued by 
Rosolowsky et al. (2003). This paper is an extension of 
our study in the Milky Way (Imara & Blitz 2010, hence- 
forth Paper I) in which we compared the velocity fields 
of GMCs to those of the local atomic gas surrounding 
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them. Our study is prompted, in part, by observations 
in external galaxies showing that high-density H I ap- 
pears to be a necesssary but not sufficient condition of 
GMC formation (e.g., Mizuno et al. 2001; Engargiola et 
al. 2003). Furthermore, it has been suggested that the 
discrepancy between observations and theoretical expec- 
tations in GMC angular momentum arises because angu- 
lar momentum has been redistributed or "shed" during 
the course of GMC evolution (Mouschovias & Paleolo- 
gou 1979; Mestel & Paris 1984; Rosolowsky et al. 2003, 
Kim et al. 2003). It seems reasonable, then, to look 
for evidence of this redistribution in the immediate envi- 
rons of GMCs. And thus, we are motivated to search for 
clues to the angular momentum problem by comparing 
the velocity fields in GMCs to those of the atomic gas 
with which they are associated. Whereas Rosolowsky et 
al. (2003) make inferences about the H I velocity ileld 
from the M33 rotation curve, in this study, we take mea- 
surements of the gradients in the atomic gas associated 
with the GMCs directly from observations. 

In the next section, we describe the observations used 
to perform this analysis. In §3, we provide the methodol- 
ogy and discuss the properties of atomic gas surrounding 
molecular clouds in M33. Implications for GMC forma- 
tion and evolution are discussed in §4, and the results 
are summarized in §5. 

2. DATA 

The data for the M33 molecular clouds are obtained 
from the Rosolowsky et al. (2003) catalog. Their work 
was a high-resolution continuation to the Engargiola et 
al. (2003) survey of 148 GMCs. Rosolowsky et al. ob- 
served 17 fields of the highest-mass clouds from their 
previous study in ^•^CO(J = — 1) using the C configu- 
ration of the BIMA array (Welch et al. 1996) in the fall 
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of 2000 and spring of 2001. The - 6" resolution of the 
observations, corresponding to a linear resolution of 20 
pc, was sufficiently high to resolve most GMCs and mea- 
sure their sizes, which was done using the prescription 
of Rosolowsky & Leroy (2006). The high spatial reso- 
lution, in combination with the velocity resolution of 2 
km s^^, also allowed them to determine velocity gradi- 
ents as small as ^ 0.01 km s~^ pc~^ in the GMCs. The 
magnitudes and directions of the gradients measured by 
Rosolowsky et al. are listed in Table 1. 

Observations of the 21-cm line in M33 were taken from 
VLA archives and reduced by Schruba (A. Schruba & A. 
LeRoy, priv. comm.). The data, taken by Thilker et al. 
in 1997, where first described by Thilker et al. (2002). 
We converted the data cube from units of flux density 
(Sv) to brightness temperature (Tb) via the Rayleigh- 
Jeans relation, 

n = 1360—^^ K = 6.055 x 10^- ^ K, 

(1) 

where 5,^ is in Jy, A = 21.1 cm is the observing wave- 
length, and bmaj x &mm = 8.3" X 7.5" is the synthe- 
sized beam size in arcseconds. The final map has 2" 
pixels. The LSR velocity range covered was —321 to 
—65 km s^^ with 1.3 km s^^ channel spacing, which was 
sufficient to include all of the gas in the disk of M33. 

3. H I PROPERTIES 

In order to compare the velocity gradients of molecu- 
lar clouds with the gradients of the material from which 
they formed, we develop a set of criteria for selecting re- 
gions of atomic gas that are associated with the GMCs 
in our sample. The challenge is to determine physically 
meaningful boundaries of atomic gas surrounding a given 
molecular cloud. Thus, for each GMC, we begin by iso- 
lating subcubes in the VLA 21-cm data set that coincide 
with the GMCs kinematically and spatially. Properties 
of the molecular clouds and the regions surrounding them 
are listed in Table 1. We use a distance to M33 of 850 
kpc for this analysis (Kennicutt 1998), and we assume a 
galactic inclination of 51° (Corbelli & Salucci 2000) when 
calculating the gas surface densities. 

3.1. Selecting H I 

Molecular clouds have well-defined boundaries at 
which a phase transition occurs between the molecu- 
lar and atomic gas and, thus, they generally have well- 
defined dimensions and masses (e.g., Kutner et al. 1977; 
Blitz & Thaddeus 1980; Blitz 1993). The neutral gas 
associated with GMCs, however, does not have such dis- 
tinct edges but rather has slower, more gradual tran- 
sitions from high- to low-surface density material. In 
trying to determine appropriate boundaries for H I asso- 
ciated with GMCs, therefore, our challenge is to define 
regions that are large enough so that most of the H I be- 
longing to a given system is included, and yet tlie regions 
must be centered close enough to the GMCs so that the 
associations are physically meaningful and we do not in- 
clude too miich unrelated material. 

There are a number of properties of H I and GMCs 
observed both in the Milky Way and M33 that help us 
develop our selection criteria. First, an H I peak is typi- 
cally found to have a velocity distribution which may be 



approximated by a Gaussian and whose line center falls 
near the CO line center of the affiliated molecular cloud 
(Andersson et al. 1991; Engargiola et al. 2003). We as- 
sume that the atomic gas associated with a given GMC 
will be at approximately the same velocity as the GMC, 
within the limits set by the velocity dispersion inherent 
to the atomic gas. Since the velocity dispersion of the 
H I gas is several times larger than that of the CO, we 
select velocities falling within ±25 km s~-^ of the CO line 
center, so that we include as much emission as possible 
from the atomic gas associated with a given GMC. We 
decide to select from a wider range of velocities than the 
previous study (i.e.. Chapter 2) because the work of En- 
gargiola et al. (2003) shows that the H I associated with 
GMCs in M33 has emission extending as much as ±35 
km from the line center of the CO emission. And be- 
cause M33 has a moderate inclination (~ 51°; Corbelli & 
Salucci), the problem of source confusion is minimized, 
allowing us to select a more generous range of velocities. 
Changing the range of velocities by a 5 - 10 km s~^ does 
not significantly change the results. This is because the 
velocity maps we create are weighted by the intensity 
(sec §3.3), so low-level emission occurring near the wings 
of the spectra docs not make significant contributions to 
the velocity maps nor, therefore, to the measurement of 
the velocity gradients. 

Identifying appropriate spatial boundaries is less clear- 
cut. Figure 1 shows that in M33, there is a spatial corre- 
spondence between GMCs and over-dense H I filaments. 
The radial accumulation length, Ra, can give us an idea 
of the extent of the atomic gas associated with GMCs. 
The size of an area which collapses to form a molecular 
cloud, Ra, can be estimated by requiring that the total 
H I mass in the region be equal to the observed mass of 
the cloud. The minimum value of Ra is determined by 
assuming that the initial configuration of atomic gas is 
a cylinder stretching to infinity in the direction perpen- 
dicular to the galactic disk: 

nRlUm = Mgmc, (2) 

where Shi is the average local surface density of atomic 
hydrogen, measured from the 21-cm data. We estimate 
Ehi for different sized regions centered on the GMCs. 
The mean value of Shi spans a narrow range of values 
for regions having a size of < 200 pc. Within 150 pc of 
molecular clouds, Shi is 8.8 Mq pc~^; within 50 pc the 
average surface density is 11.5 Mq pc~^. Figure 2 plots 
the mean H I surface density within 70 pc of GMCs ver- 
sus GMC mass. The figure shows that as GMC mass in- 
creases, the surface density of associated H I (and there- 
fore HI mass) slowly increases as Ehi oc Mq^q. We 
find that the average H I surface density in the vicinity 
of GMCs is roughly 10 AIq pc^^. This is the satura- 
tion level of atomic gas, observed in other spiral galaxies 
including the Milky Way, above which gas becomes pri- 
marily molecular (Martin & Kennicutt 2001; Wong & 
Blitz 2002; Bigiel et al. 2008). We calculate Ra for the 
slightly different estimates of Shi and find a mean value 
of ^ 70 pc. Thus, we define H I regions in the 21-cm 
data cube as 140 pc x 140 pc x 50 km s^^ connected 
regions that are centered on the locations of GMCs in 
position and velocity space. Note that the area defined 
by the accumulation radius of a GMC contains gas which 
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has filled in the space previously occupied by that which 
formed the GMC. 

We briefly note that defining H I associated with 
GMCs by identifying contiguous regions in position- 
velocity space above a given temperature threshold 
would not have been appropriate for our purposes. Given 
that the H I has such a diffuse distribution compared to 
the molecular gas bound in GMCs, it is unclear what 
such a threshold should be. Even if we were to define 
H I clouds using this method, we would still be left with 
the problem of deciding which clouds in the vicinity of 
a given GMC were associated with the said GMC and 
which which were not. Based on this rationale, we find 
it much more suitable to identify H I regions using the 
accumulation radius analysis described above. 

The top-right panels of Figures 3-8 show surface den- 
sity maps of the H I surrounding six of the 45 GMCs in 
this study. The bottom-right plot displays the average 
'^'^CO and H I spectra toward the clouds. (The left-hand 
panels are described in §3.3.) Based on the spatial and 
kinematic proximity of GMCs to local maxima in the 
atomic gas, we observe three rough "classes" of clouds. 
Most GMCs in the catalog are spatially and kinemati- 
cally coincident with local H I maxima. Clouds 1, 4, 39, 
and 45 arc representative examples — at least 1/3 of the 
area of the molecular cloud coincides with a local peak in 
the atomic gas. Twenty-nine out of 45 the GMCs, 64%, 
fall into this class. The average GMC mass in this class 
is 2 X 10^ Mq, with a dispersion of 80%. Another group 
of clouds is kinematically but not quite spatially coin- 
cident with H I maxima. GMCs in this class, including 
Clouds 16, 24 (Figure 6), 35, and 44 are located near the 
edge of a bright H I peak or sit on a filament between 
two peaks. The thirteen (29%) molecular clouds belong- 
ing to this class tend to have their mean LSR velocities 
offset from the mean H I velocity by a few km s~^. The 
typical GMC mass in this second class is nearly 3 times 
smaller than in the first: 7 x 10*^ Mq, with a dispersion 
of 30%. 

At least two GMCs in the catalog. Clouds 13 (Figure 

5) and 38, are not associated with any local H I maxi- 
mum. For both of these clouds, the center of the GMC 
is roughly 60 pc from the nearest peak in H I surface 
density, and the central velocity of the molecular and 
atomic clouds differ by at least 10 km s^^ At 2.05 x 10^ 
Mq, the total mass (molecular and atomic) in the region 
surrounding Cloud 13 is among the lowest in this study; 
only four other regions have less total mass, including 
Clouds 7, 8, 9, and 43. Cloud 38 does not appear to 
have any other characteristics, such as mass or size, that 
distinguishes it from others in the catalog. The one other 
case in which the there are no H I peaks in the vicinity of 
the molecular cknid may be Cloud 43. It has the lowest 
mass of the sample, and it is one of three GMCs farthest 
from the galactic center. However, because of its anoma- 
lously high LSR velocity of —74 km s^^, compared to a 
mean galactic value of —111 km s^^ at that location, this 
molecular cloud — which was not observed in the original 
Engargiola et al. (2003) catalog may be a false detec- 
tion. In their study of star formation in the outer parts 
of M33, Cardan et al. (2007) also found GMCs that 
were not associated with strong H I emission. Neither of 
these two clouds, located at (0l'^34'"25% 30°54'50") and 



(0l'^34'"16.9^ 30°59'31.4" ), were detected by Engargi- 
ola et al. (2003), whose study was confined to the central 
region of M33. 

3.2. H I Without GMCs 

Since over-dense regions of atomic hydrogen are gener- 
ally necessary but not sufficient for molecular cloud for- 
mation, we would like to know what the distinguishing 
features are between atomic gas harboring GMCs and 
regions of dense H I in which GMCs have not been de- 
tected. One way to address this issue is to search for ve- 
locity gradients along "non-GMC" H I filaments in M33 
and then compare the properties to regions of H 1 con- 
taining GMCs. We would like to know, for example, if we 
can find regions of H I having low angular momentum — 
that is, angular momentum comparable to that which is 
observed in GMCs? 

We select a population of 45 over-dense H 1 regions — 
corresponding to the number of GMCs in the Rosolowsky 
et al. (2003) catalog— in the M33 map along H I fil- 
aments in which GMCs have not been observed. The 
regions, displayed in Figure 9, must have a mean sur- 
face density of at least 10 Mq pc~^ within 140 pc x 140 
pc, the same size as the regions of H I associated with 
GMCs. Examples of the non-GMC surface density maps 
and H I spectra are shown in the right-hand panels of 
Figures 10 12. 

To choose a relevant range of velocities at each loca- 
tion, we begin by estimating the mean velocity at which 
the emission peaks. The mean H I velocity, vq, at a given 
location is derived from the intensity-weighted first mo- 
ment of the entire velocity distribution. For the H I as- 
sociated with GMCs, we select velocities falling within 
±25 km s^^ of the CO line center. For the non-GMC 
H I, the surface density and velocity maps arc created by 
integrating over velocities within ±25 km s^^ of vq. Ta- 
ble 2 lists the coordinates, central velocities, and other 
properties of the non-GMC H I regions. 

3.3. Velocity Gradients 

We create intensity-weighted first moment maps both 
of the 140 pc x 140 pc x 50 km s^-'^ H I regions sur- 
rounding the molecular clouds and of the atomic gas not 
observed to contain molecular clouds. Following Imara & 
Blitz (2010), we fit planes to the resulting velocity maps 
and solve for the coefficients of the fit, which define the 
velocity gradient, JIhi, and the gradient direction, Om. 
The maps of regions containing GMCs are shown in the 
upper-left panels of Figures 3-8 and, of non-GMC re- 
gions, in Figures 10 - 12. The results of the fits are given 
in Tables 1 and 2. 

We find that (1) approximately 41 out of 45 of the 
H I patches harboring GMCs have significant linear ve- 
locity gradients; the mean for this group is 0.050 ± 0.004 
km s^^ pc^^ (Figure 13); (2) the mean of the gradients in 
the non-GMC H I, 0.033 ±0.003 km s"^ pc-\ are gener- 
ally smaller than those measured in atomic gas observed 
to contain GMCs; (3) the position angles of the GMC 
rotation axes are generally not aligned with the axes of 
the surrounding atomic gas, as seen in Figure 15; and 
(4) about 62% of the H I regions associated GMCs have 
position angles that differ from the axis of the galaxy by 
more than 90° (Figure 16). 
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3.3.1. Gradient Magnitudes 

The upper-left panels of Figures 3-8 show the first 
moment maps of the atomic gas within 140 pc x 140 
pc of six of the GMCs in this study. ^ The position- 
velocity plots on the bottom-left indicate whether planes 
are good fits to the first moment maps. These plots show 
the mean velocity at a given location in the velocity map 
as a function of perpendicular offset from the rotation 
axis. When a position- velocity plot can be fitted with a 
non-zero slope, m, this indicates that a plane is good fit 
to the first-moment map and that the velocity gradient 
can be assumed to be linear. The slope of the position- 
velocity plot and the uncertainty in the least-squares fit 
of the slope, <Tm, are calculated for each gradient individ- 
ually and arc indicated in the figures; the ratio of these 
parameters is listed in Tables 1 and 2. 

Based on examination of the first moment maps and 
the position- velocity plots, we set the following crite- 
ria for degree of linearity in the gradients. Regions of 
H I having position- velocity plots < 2 am have non- 
linear, random velocity fields. Regions with m « 2-3 
cTm have "marginally" linear velocity gradients; Figiirc 7 
is an example of such a borderline case. Regions with 
TO > 3 am have unambiguously linear gradients. We find 
that 39 out of 45 H I regions associated with GMCs fall 
into this last category. Applying the same criteria to the 
non-GMC atomic gas, 43 regions have m > 3 a.^- Exam- 
ples of the velocity fields in these regions are displayed in 
Figures 10 - 12. We do not observe a significant correla- 
tion between the gradient magnitudes of molecular and 
atomic gas (Figure 14). Tables 1 and 2 list the minimum 
linear extent of the gradients, as well as their magnitudes 
and directions. 

Figure 13 displays the distribution of gradient mag- 
nitudes in the molecular and atomic gas. The mean of 
the magnitudes in the atomic gas surrounding GMCs, 
0.050 ±0.004 km s~^ pc~^, are comparable to those mea- 
sured in clouds in the Milky Way (Imara & Blitz 2010). 
On average, the magnitudes in the H I around GMCs 
are less than the typical 0.07 km s~^ pc~^ observed in 
GMCs. The average gradient observed in non-GMC 
H I is 0.033 ± 0.003 km s'^ pc'^. 

3.3.2. Gradient Directions 

Many authors have argued that the linear velocity gra- 
dients observed in molecular clouds are caused by rota- 
tion (e.g., Kutner et al. 1977; Blitz 1993; Phillips 1999). 
If this is the case, the velocity gradient magnitude of a 
cloud is a measure of its angular speed and the gradi- 
ent direction, 9, is perpendicular to its spin axis. Figure 
15 demonstrates that 53% of GMCs have position angles 
that differ from the those measured in the surrounding 
atomic gas by more than 90°. Thus, if the molecular 
clouds are rotating, slightly more than half of the clouds 
are counter-rotating with respect to the atomic gas with 
which they are surrounded. 

Rosolowsky et al. (2003) show that for a galactic kine- 
matic position angle of (pg^i = 21°, GMCs are prefer- 
entially aligned with the galaxy and approximately 40° 
are retrograde rotators. We find that localized regions 
of H I — both with GMCs and not observed to contain 

^ See the online appendix for figures of the remaining H I regions. 



GMCs have gradient directions that are generally un- 
aligned with respect to the galactic axis. Figure 16 shows 
that 62% of the H I regions with GMCs have gradient 
position angles that differ by more than 90° from the 
sense of galactic rotation. The median difference between 
H I position angles and (pgai is —90° for H I containing 
GMCs and -60° for non-GMC H I. 

Rosolowsky et al. (2003) demonstrate that GMCs ly- 
ing close together tend to have their gradients aligned. 
Similarly, we find that neighboring regions of H I contain- 
ing GMCs have comparable gradient directions, within 
1° to 15° of each other. This is not surprising, since for 
neigboring GMCs we are taking measurements in over- 
lapping regions of atomic gas. In fact, sixteen of the 
molecular clouds in this sample have accumulation radii 
which overlap by at least {1/2)Ra, raising the question 
of whether the distribution of ^hi — 4'gai, shown in Fig- 
ure 16, is biased due to over-counting. We consider this 
by taking the average (/)hi for regions with accumiilation 
radii overlapping by at least 1/2. We count such regions 
once and find that the distribution of (pm — cpgaj is not 
significantly altered. The median difference between be- 
tween (^Hi and (f)ga\, in this case, is still —90°. 

Figure 17 displays the velocity vector field of both 
groups of atomic gas. The arrows point in the direc- 
tion of increasing velocity and have sizes proportional to 
the gradient magnitudes. Along the H I filament in the 
the North-East of M33, the non-GMC H I regions tend to 
have gradient directions pointing toward the South and 
South- West. And a string of non-GMC H I regions along 
the prominent filament in the South- West Clouds 23A, 
21A, 19A, 13A, 12A, 17A, and 15A— ah have gradients 
pointing within 40° of due South. Gradient alignment 
is not always observed along filaments, however. The 
non-GMC H I regions in the nearly vertical filament at a 
right ascension of of ^ 01^33™30^ have a seemingly ran- 
dom distribution of gradient position angles. We also ob- 
serve adjacent regions of H I associated with GMCs hav- 
ing large differences in gradient directions — for instance, 
Clouds 20 and 23 located near ~ 01^34™, 30°39 have 
gradients pointing in nearly opposite directions. Over- 
all, neither population of H I appears to have a global, 
systematic pattern. 

4. IMPLICATIONS FOR CMC EVOLUTION 

Assuming that molecular clouds are rotating, measure- 
ments of their angular momentum may reveal informa- 
tion about their past. The angular momentum per unit 
mass, j, is a useful quantity for comparing the angular 
momenta in different regions having comparable mass. 
For a rotating body with a power-law density profile, the 
specific angular momentum is 

j=/30i?2, (3) 

where /3 is a constant determined by the mass distribu- 
tion, and R is the cloud radius. The constant ,6 ranges 
from 0.33 for oblong bodies to 0.5 for disks. We adopt 
the intermediate value, (3 = 0.4 (used for spherical struc- 
tures having constant surface density), for GMCs. Note 
that because we do not know the inclination, i, of a given 
cloud along our line of sight, our measurement of O un- 
derestimates the true value, Q/ sini and, hence, our mea- 
surement of j is also underestimated. The position angle 



5 



of the gradient is the direction of the total angular mo- 
mentum vector. 

Equation 3 may also be used to estimate the specific 
angular momentum initially imparted to a GMC by the 
ISM from which it forms. This depends on the process 
of GMC formation, and here, following Rosolowsky et al. 
(2003) and Paper I, wc assume that GMCs collapse or 
condense from the surrounding atomic gas via a "top- 
down" mechanism, such as a large-scale gravitational in- 
stability. In this case, Q is the velocity gradient in mea- 
sured in the galactic disk at the location of the GMC, 
the accumulation radius Ra (Equation 2) is the size of 
the region that collapsed to form the GMC, and /? = 0.5. 

Figure 18 displays the distribution of specific angular 
momenta for GMCs and the predicted j for the two pop- 
ulations of HI regions. The histograms show that, if 
the gradients observed in the molecular clouds are due 
to rotation, and assuming simple top-down formation, 
the specific angular momentum of the molecular clouds 
is less than that of the atomic material from which they 
presumably formed. Rosolowsky et al. (2003) found that 
the angular momentum expected from simple formation 
theories, including the Toomre and Parker instabilities, 
was higher than the observed value by an average fac- 
tor of 5 and by as much as an order of magnitude. Wc 
measure an average factor of 27, a median of 13, with 11 
clouds having a ratio of more than an order of magnitude. 
Figure 19 shows that jui > jcMC is always the case for 
the 36 resolved GMCs in the Rosolowsky et al. (2003) 
catalog. The reason that we are measuring larger ratios 
between the expected and observed j is because wc have 
taken direct measurements of the H I velocity gradients, 
Ctui, whereas Rosolowsky et al. used a galactic rotation 
curve (Corbelli & Schneider 1997) to estimate ilni. In 
other words, we have measured higher velocity gradients 
in the atomic gas surrounding GMCs than those deter- 
mined from the galactic rotation curve. Overplotted in 
Figure 19 are data points from our study of five molecular 
clouds in the Milky Way (Paper I). The Galactic GMCs, 
including Perseus, Orion A, NGC 2264, MonR2, and the 
Rosette, follow the trend of the M33 clouds in that the 
ratio of the specific angular momenta, jHi/iGMCi ^il- 
ways greater than unity and decreases with increasing 
icMC as (jHi/jGMc) oc Jgmc This proportionality 
is driven by the equality between jgmc and its inverse, 
as we do not observe a significant correlation between jui 
and jGMC- 

If GMCs inherit their rotation from the ambient 
H I from which they condense, one might expect their 
gradients to be much larger than observed. By conser- 
vation of angular momentum, in the absence of external 
forces, the ratio of the velocity gradient magnitudes be- 
fore and after contraction is proportional to the radii of 

the rotating bodies: fiafter/^^before = (^before/-Rafter)^- 

Thus, to conserve angular momentum, a cloud that con- 
tracts from an accumulation length of 70 pc to a typical 
GMC radius of 25 pc would have to spin up by a fac- 
tor of ^ 8. Figures 13(a), 13(b), and 14 show, however, 
that the magnitudes of GMC gradients are rarely much 
greater than those of the H I with which they are as- 
sociated. From Figure 14 and Table 1 we see that the 
average value of ilcMc/^Hi is 1.8 and reaches a maxi- 
mum of 4.6. From conservation of angular momentum, 



wc might also expect to find molecular clouds rotating in 
the same sense as the surrounding ISM from which they 
formed. As previously shown in Figure 15, however, the 
position angles in the GMCs and associated H I differ by 
at least 90° in over half of the regions. 

These key results may imply one of the following: (1) 
Since GMCs formed, external torques, such as magnetic 
fields, may have caused the redistribution of angular mo- 
mentum in the molecular clouds; (2) the H I with which 
GMCs are presently associated is unrelated to and un- 
representative of the atomic gas which originally formed 
the molecular clouds; or (3) the H I and GMCs are to 
some extent associated, but the origin of the GMC ve- 
locity gradients is not rotation. We consider each of these 
in turn below. 

1. First, given that the interstellar media of galax- 
ies are magnetized — the mean field strength of M33 is 
~ 6 ± 2 (Beck 2000) — magnetic fields potentially 
play dominant roles in the formation and evolution of 
GMCs. For instance, many authors have advanced the 
Parker instability as a mechanism for GMC formation 
(e.g., Mouschovias et al. 1974; Blitz & Shu 1980; Shi- 
bata & Matsumoto 1991). Parker (1966) showed that in 
a gravitational field, a vertically layered gas coupled to a 
magnetic field is unstable to long-wavelength perturba- 
tions. The attractiveness of this magnctoliydrodynamic 
(MHD) instability is that it has time and length scales 
that are comparable to the estimates of GMC lifetimes 
and to the distance observed between GMCs. Thus, if 
MHD effects play a significant role in the formation of 
GMCs, presumably they might also be responsible for 
redistributing angiilar momentum during later stages of 
GMC evolution. 

While Rosolowsky et al. (2003) ruled out the Parker in- 
stability as the dominant mechanism for GMC formation 
because it over-predicts the amount of angular momen- 
tum observed in GMCs, they suggested that magnetic 
braking could slow down the rotation of GMCs, since the 
Alfven speed in the ISM of M33 (~ 6 km s~^) is compa- 
rable to the timescales of certain instabilities they eval- 
uated. Yet because only 10%: of the GMCs in their cat- 
alog have values of j approaching values expected from 
j(HI), they propose that magnetic braking must occur 
early on in the lifetime of GMCs or else occurs in the 
atomic gas surrounding GMCs. Mouschovias & Paleolo- 
gou (1979) predicted a set of observational consequences 
for a rotating cloud collapsing out of an ISM threaded 
by a magnetic field that is initially perpendicular to the 
cloud's axis of symmetry. As the cloud collapses, one ex- 
pectation is that the surrounding medium in the vicin- 
ity of the cloud will be set into corotation with cloud, 
a phenomenon that we do not observe in our results. 
Mouschovias & Palcologou also predicted that magnetic 
braking is highly efficient at slowing down rotation. Ac- 
cording to their simulations, a rotating cloud with mass 
10^ Mq and density 10^ cm will have its angiilar mo- 
mentum reduced by at least 99% within 1 Myr, and the 
efficiency of magnetic braking will increase as the cloud 
contracts and as the density contrast between it and the 
surrounding medium increases. It would be worthwhile 
to find out what updated models predict for the rate of 
angular momentum loss in clouds having the masses and 
densities observed in GMCs. 

2. If the atomic gas presently surrounding GMCs is 
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unrelated to the past formation history of GMCs, the 
comparisons we are making may be invahcl. In order to 
resolve this issue we need some sort of control field with 
which to compare the H I currently harboring GMCs. It 
is for precisely this reason that we performed the anal- 
ysis, outlined in §3.2, on "non-GMC" H I regions. We 
found that the two populations of H I have distinctly 
different distributions of gradient magnitudes and posi- 
tion angles. The atomic gas associated with GMCs has a 
typical of 0.05 km s~^ pc~^ and the non-GMC H I has 
a typical value of 0.03 km s~^ pc~^. Also, the median 
difference between the position angle of the H I and that 
of the galaxy is —90° for the former population and —60° 
for the latter. These differences in fl and (j) between the 
two groups indicate that something unique has occurred 
in atomic gas associated with GMCs — presumably, some- 
thing that has either caused the formation of the GMCs 
or has resulted from their formation. 

We have already seen that the mean gradient mag- 
nitude in atomic gas associated with GMCs is greater 
than the typical gradient in non-GMC H I (§3.3). This 
implies that H I containing GMCs has higher angular 
momentum (and/or shear) than non-GMC H I. We can 
evaluate how significant the difference is between the two 
populations by investigating the parameter /3rot) the ra- 
tio of a cloud's rotational energy to its self-gravitational 
energy (Goodman et al. 1993): 

^'°'-iGM/R- 

Figure 20 shows that non-GMC H I has a narrower dis- 
tribution of /^lot and a lower mean, /3rot = 0.10 ± 0.02, 
than does atomic gas containing GMCs. Thirty-two of 
45 of the non-GMC H I regions have /3,-ot < 0.1, whereas 
only 14 of the H I regions containing GMCs have val- 
ues of /3rot < 0.1. Note that for H I regions containing 
GMCs, we calculate /3rot excluding the gravitational po- 
tential energy due to the GMCs within them (Figure 20 
[a]). This group has an average /3rot of 0.39±0.09, with 4 
clouds having /^rot > 1 (Clouds 1, 2, 17, and 43), imply- 
ing that these regions may have a considerable amount 
of rotational energy. Note that for H I regions contain- 
ing GMCs, we calculate /3rot excluding the gravitational 
potential energy due to the GMC within it (Figure 20 
[a]). If we calculate ftot including all of the gas within a 
given region, the average for H I containing GMCs goes 
down to 0.24 ± 0.04. The regions listed above still ap- 
pear to have a considerable amount of rotational energy 
in comparison to their gravitational energy: /Jj-ot = 0.90, 
0.75, 0.88, and 1.1 for clouds 1, 2, 17, and 43, respec- 
tively. Clouds having large values of j^rot can potentially 
become stable against gravitational instabilities. It is 
difficult to sec how GMCs coiild have formed under siich 
conditions. These observations suggest that between the 
time prior to the onset of GMC formation to the time af- 
ter formation, processes occur that increase fi observed 
in the H I associated with GMCs. An alternative is that 
GMCs preferentially form from gas having high fi, possi- 
bly in regions that are unstable to gravitational collapse. 

It is also worth pointing out that H I regions falling into 
Classes 2 and 3 (see §3.1) regions where the GMCs are 
not quite spatially coincident with local H I maxima — 
tend to have lower gradient magnitudes (f2 ~ 0.04 



km s~^ pc^^ with a dispersion of 40%) than regions in 
Class 1 (0.06 km s^^ pc^^ with a dispersion of 50%). 
We apply a t-test and find that the difference between 
the two classes is significant to the 95% confidence level. 
However, molecular clouds belonging to cither class do 
not appear to have different distributions of il; that is 
GMCs in both Class 1 and Class 2 have the same aver- 
age gradient (~ 0.07 km s~^ pc~^). This suggests that 
activity within the GMCs may have more of an impact 
on the surroundings than the surrounding environment 
has on the velocity fields of the GMCs. 

3. Given that the gradient directions of the GMCs 
and associated atomic gas are generally not aligned, this 
raises the possibility that the linear velocity gradients 
observed in molecular clouds may not be caused by rota- 
tion. Burkert & Bodenheimer (2000) demonstrated that 
turbulent velocity fields can also cause linear gradients. 
They found that the gradient magnitude of turbulent 
cores scales with size as 51 cx R^^ . As shown in Fig- 
ure 21, fitting a power-law relationship to the gradients 
observed in the H I surrounding GMCs as a function of 
the accumulation radius, we find fim oc ^ The 

relationship between Oqmc a-nd size for the 36 resolved 
GMCs is Ogmc oc i?-"-3±"-2. Figure 21 also displays a 
combined fit for atomic gas and GMCs, which is in good 
agreement with the Burkert & Bodenheimer (2000) re- 
sult. 

In our previous study (Paper I), we found that star 
formation activity provides a reasonable explanation of 
the gradients observed in a small sample of Milky Way 
molecular clouds. In at least three out of five cases, in- 
cluding Orion A, NGC 2264, and the Rosette, the veloc- 
ity fields may have been produced by expansion driven 
by highly energetic stellar winds. The location of H II re- 
gions corresponds to the highest velocities of the molec- 
ular material in these clouds. Though the sample was 
too small to make a firm conclusion, it appears that the 
H II regions within these molecular clouds may also be 
influencing the morphology of the associated atomic gas. 
In NGC 2264 and the Rosette, for instance, there are 
local peaks of H I that appear to have been swept up by 
expanding molecular gas. Following up this study of M33 
with an analysis comparing the locations of H II regions 
to the kinematic properties of GMCs could provide useful 
insight to how the velocity fields of the latter evolve. 

5. SUMMARY 

We have presented a detailed comparison of the veloc- 
ity fields in 45 GMCs detected in M33 by Rosolowsky 
et al. (2003) using BIMA observations and the atomic 
hydrogen with which they are associated. Using high- 
resolution VLA 21-cm observations to create surface den- 
sity and intensity-weighted first-moment maps of the H I, 
we also compared the properties of atomic gas contain- 
ing molecular clouds with atomic gas in which molecular 
clouds have not been detected. Based on our measure- 
ments, including the velocity gradient magnitudes and 
directions of these regions, we make the following con- 
clusions: 

1. The average surface density of atomic hydrogen as- 
sociated with the GMCs is ~ 10 Mq pc~^, similar to the 
saturation level of H I above which gas becomes primar- 
ily molecular in other galaxies, including the Milky Way 
(Martin & Kennicutt 2001; Wong & Blitz 2002; Bigiel et 
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al. 2008). A power-law relationship exists between the 
H I surface density and GMC mass: Shi oc M°-^'^(f 

2. We observe three categories of GMCs, based on their 
proximity to local peaks in the atomic gas. The major- 
ity of GMCs (64%) coincide spatially and kincmatically 
with local H I peaks. Twenty- nine percent of GMCs arc 
located near the edge of an H I peak, or sit between two 
peaks. Clouds in this category tend to have their mean 
LSR velocities offset from the mean H I velocity by a 
few km s~^. The remaining three clouds (7%) are not 
associated with H I maxima and the mean velocities of 
the molecular cloud and atomic gas are offset by ~ 10 
km s~^. 

3. Thirty-nine out of 45 of the H I regions in 
the vicinity of GMCs have linear velocity gradients of 
around 0.050 km s~^ pc~^ and spanning 0.013 to 0.13 
km s"'^ pc"-'^. If GMCs arc rotating and initially inher- 
ited their angular velocity from the surrounding H I, con- 
servation of angular momentum would require that the 
angular speed of a typical GMC {R « 25 pc) speed up by 
a factor of 8. The average value of Ogmc/^hi, however, 
is only 1.8. Magnetic braking has been used by some au- 
thors to explain how the slowing of GMC rotation may 
have occurred. But if magnetic braking is as efficient at 
slowing rotation as predicted (e.g., Mouschovias & Pale- 
ologou 1979), we would expect most GMCs to have much 
lower values of Ogmc — as much as an order of magni- 
tude less — than observed. Alternatively, magnetic brak- 
ing may be less efficient, or it may operate over longer 
timescales. 

4. Fifty-three percent of the molecular clouds have 
gradients whose directions differ from the gradient direc- 
tion in the local H I by more than 90°. If the gradients 

in the GMCs were caused by rotation, this implies that 
over half of them arc counter-rotating with respect to the 
atomic gas with which they arc surrounded. This mea- 
surement is difficult to reconcile with the notion that 
gradients are due to rotation and that GMCs form in a 
simple top-down matter. 

5. Gradients in the atomic gas associated with GMCs 
generally have larger magnitudes than expected from 
galactic differential rotation alone. Also, 62% of these 
regions have gradient position angles that differ from the 
sense of galactic rotation by more than 90°. 

6. We examined the properties of high-density H I re- 
gions in which molecular clouds have not been detected 
and found that they have a lower range of velocity gra- 
dients, ^ 0.03 km s~^ pc^^, than regions where GMCs 
are observed. This suggests that something occurs dur- 
ing the course of GMC evolution that may increase the 
shear of the atomic gas. Neither population of atomic 
gas has gradient directions that are preferentially aligned 
with the kinematic position angle of the galaxy, nor did 
we find a correlation between gradient magnitude and 
direction in either population. 

7. A power-law relationship exists between gradient 
magnitude and size in both the molecular clouds and the 
H I surrounding them. For GMCs, r^GMC « i?-o.3±o.2. 
for H I, r^Hi oc i?^"'^^"'^, where Ra is the accumulation 
length. The combined relationship is 17 oc _R^°-^^"-^, con- 
sistent with what Burkert & Bodenheimer (2000) found 
for the velocity fields of turbulent molecular cores. 

8. Our analysis raises considerable doubt to the hy- 



pothesis that the origin of GMC velocity gradients is ro- 
tation. Alternative explanations worth exploring include 
turbulence, shear, and star formation activity. 
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TABLE 1 
Cloud Properties in M33 
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TABLE 1 

The cloud numbers correspond to those in the Rosolowsky et al. (2003) c;ataloc;. The gradient magnitudes and directions 

ARE fl AND 6, RESPECTIVELY. ThE DIFFERENCE BETWEEN THE H I AND GMC GRADIENT DIRECTIONS IS ^hi — ^GMC- THE EXTENT OF THE 
GRADIENT IS LISTED FOR H I REGIONS HAVING LINEAR GRADIENTS. ThE GRADIENT MEASUREMENTS OF CLOUDS IN WHICH LINEAR GRADIENTS 

HAVE NOT BEEN DETECTED {\m\/am < 3) ARE LISTED IN BRACKETS. 
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TABLE 2 
Properties of Non-GMC H I in M33 
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TABLE 2 

Properties of H I regions in which GMCs have not been detected. 
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R,A. (2000) 



Fig. 1. — M33: Grayscale image in units of Mq pc~^ of the 21-cm emission of the central 45' X 45' field. Molecular clouds are overlaid 
with area scaled to mass. Nearly all (93%) GMCs lay in regions of high-density H I. The galactic mean of Shi is roughly 4 Mq pc~^, 
while the mean value in the vicinity of GMCs is 10 Mq pc^^. 




Fig. 2. — Mean local Shi versus GMC mass. The solid line is the least-squares fit for the plotted points, which show the mean H I surface 
density within 70 pc of GMCs. The dotted and dashed lines are the least-squares fits for Shi within R < 150 pc and Shi within i? < 50 
pc, respectively. The average surface density of atomic gas surrounding GMCs is 10.2 ± 0.4 Mq pc~^. The H I surface density increases 
slowly with increasing GMC mass as Shi oc ^qmc'^ '^^' '^^^^^ the power-law exponent is the mean of the three least-squares fits and the 
uncertainty is the l-tr spread. 
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Fig. 3. — Cloud 1: The top left figure shows the intensity-weighted first moment map of the H I with the gradient axis overlaid. The 
velocity range of the map is indicated in the top left corner in units of km s^^; red represents the maximum speed. Below this figure is 
a plot of the central velocity at a given location in the first-moment map versus perpendicular offset from the gradient axis; the linearity 
of the plot indicates that a plane is a good fit to the first-moment map; the radial extent of the GMC is demarcated by the horizontal 
blue line. The top right figure is a surface density map of the H I overlaid with a circle matching to the size of the associated GMC; the 
10 M© pc~^ contour is overlaid in yellow. The range of H I surface densities displayed in the map are in the top left corner in units of 



Mq pc ^, and the total H I mass in the region is written in the bottom loft corner, 
(black) and CO emission (red) toward the region. 



Below is a plot of the average spectra of If I emission 
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Fig. 4. — Cloud 4: Same as Figure 3. 
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Fig. 5. — Cloud 13: Same as Figure 3, accept that the slope of the position-velocity plot is close to zero, indicating that a plane is not a 
good fit to the first-moment map. Cloud 13 is an example of a region that does not exhibit a linear gradient and in which the CMC does 
not coincide with a local maximum in the H I. 
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Fig. 6. — Cloud 24: Same as Figure 3. 
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Fig. 7. — Cloud 39: Same as Figure 3. 
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Fig. 8.— 



Cloud 45: Same as Figure 3. 



18 




R.A. (2000) 



Fig. 9. — M33: Grayscale image in units of AIq pc ^ of the 21-cm emission of the central 45' X 45' field. The locations of H I regions 
containing molecular clouds are overlaid in orange, and the locations of H I regions without observed molecular clouds are overlaid in blue. 
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Fig. 10. — Cloud lA: Same as Figure 3, accept for a region in which GMCs have not been observed. 
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Fig. 11.— Cloud 23A: Same as Figure 10 
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Fig. 12.— Cloud 30A: Same as Figure 10 
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Fig. 13. — Gradient magnitudes for (a) GMCs, (b) H I clouds containing GMCs, and (c) H I clouds without observed GMCs in M33. 
Clouds having a position angle differing from the galaxy by more than 90° are given negative values. The hatched portions of the histograms 
in (b) and (c) represent regions having non-linear linear gradients. 
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Fig. 14. — Gradient magnitudes in the atomic gas versus GMC gradient magnitudes. There is no significant correlation between and 




Fig. 15. — Comparison of H I and GMC position angles. // the velocity gradients indicate rotation, then most GMCs (~ 53%) are not 
rotating in the same sense as the associated H I. The hatched portions of the histograms in (b) and (c) represent regions having non-linear 
gradients. 



24 




^ I I ^ 

0! \ \ \ \ \ 

-180 -90 90 180 



Fig. 16. — Position angles of (a) GMC-harboring H I and (b) non-GMC H I, with respect to that of M33. The hatched portions of the 
histograms in (b) and (c) represent regions having non-linear gradients. 




0l'35""0= 34"30' 0' 33"'30'' 0' 32'"30= 

R.A. (2000) 



Fig. 17. — M33: The directions of the gradients in the atomic gas are plotted for H I regions containing molecular clouds (orange) 
and for H I regions without observed molecular clouds (blue). The arrows point in the direction of increasing velocity and have lengths 
proportional to the gradient magnitude. The gradient directions of the individual velocity fields where GMCs are observed (or, where 
GMCs may potentially be in the process of forming) do not appear to make up a large-scale, systematic pattern. 
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Fig. 18. — Distribution of specific angular momentum for (a) GMCs and for (b) H I clouds containing GMCs. 
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Fig. 19. — Ratio of specific angular momenta in atomic gas and GMCs, jm/iGMCi versus specific angular momentum in the 36 resolved 
GMCs. This plot shows that jni > Jgmc is always the case. The dotted line shows the least-squares fit to the data; (jm/jcMc) oc 
ioMcJ^" " ' ■ "^^^'^ median j'hi/jgmc is 13 and the average is 27. Data points for Milky Way GMCs are overplotted in filled circles (but are 
not included in the fit). 
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Fig. 20. — Distribution of the /3rot parameter for (a) H I containing GMCs and (b) non-GMC H I. 
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Fig. 21. — Gradient magnitudes observed in GMCs (red) and associated atomic gas (black) as a function of size. The data are averaged 
in bin sizes of AR = 0.1 dex. The lines indicate least-squares, power-law fits to the data. For GMCs, f2GMC R~^-^^^-^ , and for 
H I,f2Hi oc R^*^'^^"'^, where Rj^ is the accumulation radius. For GMCs and H I combined is C oc R^^'-^^'^-^ , which is the relationship 
found by Burkert & Bodenheimer (2000) for turbulent molecular cores. 
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APPENDIX 

The following figures display the velocity maps, position-velocity plots, surface density maps, and spectra of the 
H I regions associated with the remaining molecular clouds not displayed in our paper on giant molecular cloud 
(GMC) angular momentum in M33 (Imara, Bigiel & Blitz 2010). If for a given region the GMC was resolved in the 
Rosolosky et al. (2003) catalog, a circle matching to the size of the GMC is overlaid on the surface density map. 
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Fig. 22. — The top left figure shows the intensity-weighted first moment map of the H I with the gradient axis overlaid. The velocity 
range of the map is indicated in the top left corner in units of km s~^; red represents the maximum speed. Below this figure is a plot of 
the central velocity at a given location in the first-moment map versus perpendicular offset from the gradient axis; the linearity of the plot 
indicates that a plane is a good fit to the first-moment map; the radial extent of the GMC is demarcated by the horizontal blue line. The 
top right figure is a surface density map of the H I overlaid with the 10 Mq pc~^ contour in yellow. The range of H I surface densities 
displayed in the map are in the top left corner in units of Mq pc~^, and the total H I mass in the region is written in the bottom left 
corner. Below is a plot of the average spectra of HI emission (black) and CO emission (red) toward the region. 
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Fig. 23.— See Figure 22. 
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Fig. 24.— See Figure 22. 
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Fig. 25.— See Figure 22. 
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Fig. 26.— See Figure 22. 



32 



Velocity Map 



-162:-150 



-110 



Cloud 8 




40 35 30 25 20 
RA (offset in arcsec) 



Surface Density Map 




40 35 30 25 
RA (offsett in arcsec) 



-145 



-150 



I -155: 



I -160- 



-165 - 



-170 



dv/ds = -7.99+ 1.41 




-100 -50 50 

Offset from axis [pc] 



100 




-190 -180 -170 -160 -150 -140 -130 
Velocity [km/s] 



Fig. 27.— See Figure 22. 
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Fig. 28.— See Figure 22. 
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Fig. 29.— See Figure 22. 
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Fig. 30.— See Figure 22. 
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Fig. 31.— See Figure 22. 
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Fig. 32.— See Figure 22. 
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Fig. 33.— See Figure 22. 
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Fig. 34.— See Figure 22. 
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Fig. 35.— See Figure 22. 
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Fig. 36.— See Figure 22. 
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Fig. 37.— See Figure 22. 
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Fig. 38.— See Figure 22. 
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Fig. 39.— See Figure 22. 



45 



^ 140 



o 
Q 



Cloud 22 

Velocity Map Surface Density Map 





355 350 345 340 335 330 
RA (offset in arcsec) 



355 350 345 340 335 330 
RA (offsett in arcsec) 



-200^ 
-202 - 
-204 - 



dv/ds = -3.54 ± 0.77 



3 -206 



> 



-210 
-212 
-214 



-100 -50 50 

Offset from axis [pc] 



100 



to 

P9 1 








-240 -230 -220 -210 -200 -190 -180 
Velocity [km/s] 



Fig. 40.— See Figure 22. 
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Fig. 41.— See Figure 22. 
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Fig. 42.— See Figure 22. 
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Fig. 43.— See Figure 22. 
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Fig. 44.— See Figure 22. 
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Fig. 45.— See Figure 22. 
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Fig. 46.— See Figure 22. 
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Fig. 47.— See Figure 22. 
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Fig. 48.— See Figure 22. 
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Fig. 49.— See Figure 22. 
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Fig. 50.— See Figure 22. 
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Fig. 51.— See Figure 22. 
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Fig. 52.— See Figure 22. 
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Fig. 53.— See Figure 22. 
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Fig. 54.— See Figure 22. 
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Fig. 55.— See Figure 22. 
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Fig. 56.— See Figure 22. 
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Fig. 57.— See Figure 22. 
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Fig. 58.— See Figure 22. 
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Fig. 59.— See Figure 22. 
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Fig. 60.— See Figure 22. 



